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1. INTRODUCTION 

The defining property of a superconductor is the stiffness of its phase. In conventional low- 
T c superconductors this stiffness is so large that in order to directly observe its finite value 
one has to create a region where it is locally suppressed, the so-called superconducting 
weak link. The resulting weak link region gives rise to a variety of Josephson effects (for 
a review, see JIJ), in which the superconducting phase becomes a measurable quantity. 

In the cuprate superconductors the sensitivity of the Josephson effect to the phase has 
played a decisive role in establishing a phase- sensitive test of the (now well established) 
unconventional ci-wave symmetry of their pairing state (for a review, see ||). 

The purpose of this paper which is based on the results of our measurements of the 
current-phase relation in the cuprates p, f|, [5], |J is to argue that besides providing a 
phase-sensitive test of the pairing symmetry, a quantitative analysis of the Josephson 
effect in the cuprates provides new insights into the microscopic physics of the cuprates, 
in particular into the reduced phase stiffness of the cuprates. 

Our experiments were carried out on two types of Josephson junctions: [001] tilt 
grain boundary junctions || |, || and c-axis junctions between a cuprate and a low-T c 
superconductor ||. Let us start with a brief discussion of these two types of junctions. 

1. Grain boundary junctions have been studied since the very beginning of the cuprate 
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research (for a review see [|7|]). In particular, it has been realized very early on that grain 
boundaries are the current limiting element of polycrystalline materials and as such are 
to be avoided in large current applications. However, after having been used in the 
phase sensitive experiments [fj, the situation changed completely when it was realized 
that cuprate grain boundary junctions might form the basis of a completely new type of 
superconducting electronics, based on the use of the so-called 7r-junctions. 

In an idealized picture, a grain boundary junction can be thought of as a planar 
interface of two grains. The junction is characterized by two angles 9\ and 9 2 between the 
normal to the interface and the crystallographic axes in the grains, see Fig. [j]. 

It is an experimental fact that the junction transparency T> depends predominantly 
on the misorientation angle 9 = 9 2 — 9\. The microstructure of the grain boundary is 
quite complicated and at the moment there exists no generally accepted picture of it. 
In this paper we adopt the point of view advocated by Hilgenkamp and Mannhart |7]], 
who view the junction as a locally underdoped (and insulating) region. The width of this 
insulating barrier is supposed to be proportional to 9, therefore leading quite naturally 
to the experimentally observed scaling T>(9) = exp(—9/9 ) with 9 « 5°. In this paper 
we shall be concerned mostly with grain boundary junctions with misorientation angle 
9 = 45° for which V m e~ 9 « 10" 4 . 





0\ 






^0 





Figure 1: Schematic drawing of a [001] tilt grain boundary junction. 

We believe that grain boundary junctions with 9 > 9 are to be described as featureless 
junctions in the tunnel limit. An alternative explanation of the very small transmission 
probability T>(9) could invoke the presence of pinholes in the barrier, whose density would 
decrease exponentially with 9. We believe that measurements || |8| of the relation between 
the superconducting current I and the phase difference ip across the junctions with 9 = 24° 
and 36° have ruled out this possibility, since the results were fully consistent with the 
tunnel limit prediction 1(f) oc simp and no trace of higher harmonics indicative of small 
barrier weak links were found, see Fig. 

2. The other type of junctions whose current-phase relation we have studied are the 
so-called c-axis junctions between the cuprates and low-T c superconductors. In this type 
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of junctions the interface is parallel to the Cu02 planes and the current flows along the c 
axis of the cuprates. Because of the layered structure of the cuprates the interlayer forces 
are presumably much weaker than those within the layers, and therefore much cleaner 
and better defined interfaces are to be expected for this type of junctions, when compared 
with the grain boundary junctions. 




Figure 2: Current-phase relation at T = 15 K for symmetric GBJJs with 9 = 24° (dotted 
line), 36° (dashed line) and 45° (solid line). The inset shows the scaling j c oc exp(— 6/ 9$) 
of the critical current density with 9q = 5.6°. Taken from ||. 



2. THEORY 

2.1. GRAIN BOUNDARY JUNCTIONS. Our main interest in this paper regards the 
current-phase relation I(ip), which is on general grounds assumed to be a 2tt periodic 
function. Moreover, in absence of magnetic fields we expect that J(— </?) = —/(</?) and 
therefore we can write ]]J 

I((p) = l\ sin ip + I 2 sin 2<p + . . . . (1) 

The c?-wave symmetry of the pairing state in the cuprates requires that odd harmonics 
Ji,/ 3 ,... change sign with 9^ — > 9i + tt/2. If we keep only the lowest-order angular 
harmonics, we can therefore write || 

It = I c cos 20i cos 20 2 + I s sin 29 x sin 20 2 . (2) 

The coefficients J c , I s are functions of the barrier strength, temperature T, etc. In partic- 
ular, if the tunneling is allowed only for a narrow range of impact angles around normal 
incidence, only the I c term survives. This situation has been considered in the early paper 
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by Sigrist and Rice |TD| , who suggested to use the first term of Eq. |2| as a simple approxi- 
mate formula which respects the symmetries of the Ambegaokar-Baratoff like expression 
(at T = 0) for the critical current, 

7rA fc Ay / |A fc - Ag| \ 

: K i a i I i a i ^ c cos26'icos26'2, (3) 



eR N (\A k \ + \A q \Y 

where Rn is the resistance of the junction in the normal state, A*., A q are gap functions 
for the case of normal incidence in the two superconducting electrodes, and K is the 
complete elliptic integral. 

On the other hand, the even harmonics J 2 , J 4 , . . . are not forced by symmetry to change 
sign and therefore we neglect their weak angular dependence (except for the strong depen- 
dence on the barrier transparency T>{9)). This then implies (as first noted in |IT| ) that for 
the so-called asymmetric 45° grain boundary junctions (i.e. 9\ = 45°, 62 = 0°), the first 
harmonic I\ vanishes by symmetry and the current-phase relation can be approximated 
by I(ip) = I2 sin 2 99. 

In real life the junction interface meanders around its mean (nominal) direction. In 
that case we wish to interpret the function I\{9\, 82) as a relation between the local critical 
current density and the local interface geometry. Structural studies of grain boundaries 



12| show that the interface is typically faceted with a typical size of a facet ~ 100 nm. 
Since this length scale is much larger than both relevant electronic length scales, the Fermi 
wavelength and the coherence length, we believe that the macroscopic first harmonic can 
be estimated by a simple averaging of Eq. ^| along the interface. Making use of such a 
procedure, if the local junction geometry is described by 6i + x{ x )i h °f a rough interface 
is given by || 

h = \{h + Q cos2£ + |(/ c - I s ) cos2(^ + 9 2 ), (4) 

where x = (cos 4%) is an interface roughness parameter and (...) denotes an average 
along the interface. Note that in the completely rough limit x = 0, I\ depends only on 
the misorientation angle 9 and not on the individual angles 9% 0, since in that limit the 
notion of the direction of the interface is meaningless. 

On the other hand, since I2 is not forced by symmetry to depend strongly on 9i (except 
for the dependence through the barrier transmission on the misorientation angle 9), we 
expect that I 2 will exhibit only a weak dependence on the surface roughness x. 

Summarizing the above symmetry arguments, it follows that d-wave symmetry and 
sufficient interface roughness (small x) imply enhanced I2/I1 ratios for all 45° junctions. 

However, while very robust, the symmetry arguments do not provide estimates of the 
ratio In what follows we attempt a more quantitative analysis of the Josephson 

effect. The increased level of detail has its prize, however: we have to introduce a mi- 
croscopic model of the junction. Two such models have been discussed in the literature: 
a perfectly flat interface which leads to the presence of midgap states in the case of an 
impenetrable barrier ||13|| , and a model emphasizing the roughness of the interface and 
the subsequent development of a nontrivial pattern of currents flowing along the interface 



2.1.1. IDEALLY FLAT JUNCTIONS. Let us consider first the 'flat' scenario. The key 
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role in this approach is played by the concept of the so-called midgap states, which 
should form at surfaces of anomalous superconductors. Theory predicts that the number 
of midgap states should be maximal for (llO)-like surfaces of d-w&ve superconductors 
and no such states should form at (100) and (OOl)-like surfaces |13[. The presence (and 
the nontrivial surface orientation dependence) of such midgap states in the cuprates is 
well documented by now by STM spectroscopy |T5[ and by grain boundary tunneling 
spectroscopy fl6fl . It has been also suggested that such bound states are at the origin of 
the observed nonmonotonic temperature dependence of the penetration depth |T7[] . 

Let us turn now to the study of the influence of the midgap states on the Josephson 
effect in the cuprates, following the early proposal of Tanaka and Kashiwaya ]l8|| (for 
reviews, see |19|, |20||). Within the simplest description, we consider a circular Fermi line, 



specular reflection at the interface, and a constant order parameter in the superconducting 
grains. Because of the translational symmetry along the junction, the scattering problem 
factorizes into a set of independent problems for each given impact angle. 

The energy of the Andreev bound states is governed by the probability D(6,d) of 
transmission for impact angle The function D{9 1 , d) depends on the details of the 
barrier, such as its width and height, none of which are known reliably in the present 
context. The only direct experimental measure of the barrier transmission is the resistance 
per square of the junction in the normal state, R a = (h/2e 2 )(ird/kFT>), which provides 
us with one moment of the function D{9,&), V{6) = £ /2 d$ cos $£>(#, tf). 
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Figure 3: Trajectories of the particles for specular reflection on a flat junction interface. 
The impact angle $ is the same for both sides of the junction. 

Let us calculate the energy of the Andreev bound state corresponding to the incoming 
(outgoing) trajectories in the left and right grains denoted 1 and 2 (3 and 4), respectively 
(see Fig. |3|). To this end, let us denote the order parameter on the trajectory i as Aj = 
|Ai|5i, where Si = ±1 is the sign. The qualitative character of the bound state depends 
on the relative value of the signs Si. 

Before proceeding, let us recall that for an impenetrable barrier with D — > 0, a midgap 
state is formed in the left (right) grain for Si = —S3 (S2 = —S4), i.e. when the sign of 
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A changes under reflection at the interface [TBI. In the tunnel limit, the energy of the 



Andreev bound state can belong to one of the following three classes: 

(i) The signs of A change under reflection on both sides of the junction, Si = —S3 
and S 2 = — S4. In that case the two midgap states on both sides of the junction are split 



by the finite value of D [21 
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(5) 



Note that e oc y/D, as is usual for degenerate perturbation theory. 

(ii) The midgap state is formed only on one side of the junction, the left one for 
definiteness (Si = —S3, S 2 = £4). In that case there is only one anomalous Andreev 
bound state, |23fl : 

[AiA 3 | 



e(9,#,<p) = D(9,#)- 



rSiS 2 sin (p. 



(6) 



|Ai| + |A 3 | 

(iii) There are no midgap states on both sides of the junction, i.e. Si = S 3 and 
5*2 = S4. In this case there are no amalous Andreev bound states and the junction 
resembles qualitatively that between two s-wave superconductors. Note that this type of 
processes is always realized for 1? = 0°. 

Once the quasiparticle energies are known, the free energy due to Andreev levels per 
unit area of the junction reads 



Fn(<p) 



kpT 
2ttc/" 



E 



tt/2 
-tt/2 



d$ cos i? In 



1 + exp 



T 



(7) 



where d is the average distance between the Cu02 planes and the index n numerates all 
^-dependent energy levels for the impact angle According to general principles, the 
current-phase relation I(ip) can be determined from the ^-dependence of the total junction 
free energy, / = (27r/$o)<9-F/<9<£> and the density of supercurrent across the junction reads 
j(cp) = (2n/Q )dF n /d(p. We emphasize that it is the total ^-dependent free energy 
density Fn(ip) (and not only the contribution of the Andreev bound states F^-) which 
determines j((p). 

Before proceeding we should like to point out that, in general, a finite phase difference 
ip across the junction leads to a nonzero current along the junction [p^| . This current 
generates a magnetic field which has to be screened by the Meissner currents in the 
superconducting electrodes [P^[ , thus leading to an additional term in the free energy, 

The spontaneously generated current along the interface can be calculated as follows. 
Since all four momenta involved in interface scattering with impact angle d (see Fig. |3|) 
have the same momentum k y along the junction, the y component of the current density 
carried by the Andreev bound state with energy e(i9) is 



X 



m 



(\u(x, ^)| 2 + \v(x, 



where u(x) and v(x) are the two components of the Andreev level wavefunction. Since 
u(x),v(x) decay exponentially on the distance ~ ^ from the interface (where ^ ~ Hvf/A 
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is the coherence length), also the net current along the interface 

3 v i?) = ^E dVcos#j y (#,x)f(e n (#)) (8) 

Z7T ~ J-tt/2 

is localized within £ from the interface. From Eq. [8] it follows that the linear density of 
the total current flowing along the grain boundary per Cu0 2 layer, J = I^Ljy( x )^ * s 

J o = ^E r £*nMf(e n (4)), (9) 

tlU „ J-n/2 2/K 

where Ep is the Fermi energy. In deriving Eq. ^ we have made use of the fact that the 
Andreev level wavefunctions are normalized, dx(\u(x)\ 2 + |t>(x)| 2 ) = (Ld)^ 1 . 




Figure 4: Schematic drawings of asymmetric (a) and symmetric (b) 45° grain boundary 
junctions. Shaded regions correspond to a positive order parameter. For asymmetric 
junctions, impact angles for which Si £2 = —1 are denoted by the hatched regions. For 
symmetric junctions the hatched regions correspond to trajectories with midgap states 
on both sides of the junction. 

Asymmetric 45° junctions. In what follows, let us apply the above formalism to the 
two special cases of interest to us, namely the asymmetric and symmetric 45° junctions. 
For asymmetric junctions, the trajectories for all impact angles are of type (ii), see Fig. |J 
Since in this case Ai($) = — A 3 (?9) = A sin 2$ and A 2 (i9) = A 4 (-$) = A cos 2$, the energy 
of the Andreev levels is 

£($) = —D($) I sin 2^1 sign(sin4$) sin<£>. 

For the sake of simplicity we consider a simple model for the barrier transmission function, 
D(d) = D(0) for \&\ < $ < vr/4 and D(ir/A) < D(0) otherwise. In this case V « 
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D(0) sin^o and the typical energy of the bound states ~ VA sets the energy scale of the 
problem. In the two extreme limits T -C VA and T 3> VA, respectively, we find 



tf k F VA 
An d 
k F V 2 A 2 

48^ Td 



sm ip\ 



cos 2(p. 



(10) 



From Eq. ^ it follows that Jo depends only on the population of the Andreev levels. Let us 
for definiteness consider the case ip > 0. Depending on the impact angle the trajectories 
can be classified into two distinct groups: The hatched regions in Fig. |] correspond to 
SiS 2 = —1, whereas for the remaining trajectories SiS 2 = 1. Thus the energy of the 
Andreev levels in the hatched (not hatched) regions is negative (positive). Therefore at 
T = only the Andreev levels in the hatched regions are occupied and from Eq. |9] we 
find Jo = 0. This effect has been called super screening in the literature ||24j| . At finite 
temperatures T ^> D{tt/A)A an explicit calculation leads to a finite current along the 
junction, 



Jn 



^QekpVF VA sin <p 



4tt d 3T + VA\smp\ 

The surface currents due to the Andreev bound states generate magnetic fields which 
have to be screened by the Meissner currents in the superconducting grains. As pointed 
this leads to an increase of the total junction energy, whose magnitude we 
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now estimate in the simple case of a two dimensional interface (in the y-z plane) between 
3D superconducting grains. Of course, most experiments are carried out on thin films 
and thus the 3D geometry is not fully appropriate. However, since the film thicknesses 
are typically ~ 1000 A which is comparable to the bulk penetration depth, our estimates 
should remain qualitatively correct. 

Since the current due to Andreev bound states flows on scale £, we model its distri- 
bution as j = (0,j, 0) with j(x) = J 5(x). Within the London theory (assuming bulk 
penetration depth A) , the distribution of the Meissner current and of the magnetic field 
can be readily calculated and the results are shown in Fig. |5|. The sum of the energy of the 
magnetic field and of the kinetic energy of the superflow (per unit area of the interface) 
is easily seen to be F^f = yU JoA/4. 

Therefore the y?-dependent part of the total free energy density F n (<p) = F^{ip) + 
3 '(<p) for T « DA and T » VA reads 



f m 



A k f 



•A 



47T d 

d k F V 2 A 2 
48^T T 



—T> \ sin<£>| + 



sm (p 



A V | sin ip\ + a 



1 - 



gjj hv F /T 
12 A 



COs2<y9, 



11^ 



respectively, where £ = Hvf/A and a = 3T/VA. Several points regarding Eq. [Ll| are to 
be stressed. First, in the low temperature limit, our result agrees qualitatively with |25 
in that the free energy is minimized for <p = or ip = it, if V < #g£/A (i.e. the bound 
state energy gain is smaller than the Meissner energy loss). In the opposite limit when 
the bound state energy dominates F n , the free energy is minimized for ip = ±7r/2. 
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Figure 5: The distribution of supercurrent density (a) in units of jo = Jo/^ an d magnetic 
field (b) in units of B Q = ^lqJq for symmetric 45° grain boundary junctions. (c,d): the 
same for asymmetric 45° junctions. In all figures we have used £/A = 0.1. 



From the normal state resistivity of 45° junctions we estimate T> ~ 10~ 4 . If we take 
£ ~ 30 A and A ~ 1500 A relevant for the cuprates, we estimate that for not too strongly 
angle dependent tunneling, the Meissner energy should dominate F n at low temperatures. 
We should like to point out that by no means this implies trivial physics. Just on the 
contrary, as shown in Fig. || dominant Meissner energy leads to a doubled periodicity of 
I ((f)- This is qualitatively similar to the well studied case when dominates. However, 
there are two major qualitative differences with respect to the case of dominant bound 
state energy: (i) the (degenerate) ground state corresponds to a state when no currents 
flow along the junction, and (ii) the current-phase relation becomes steep in the minima 
(not maxima) of the junction energy, see Fig. 

At temperatures above the bound state energy T>A (which is presumably the case 
down to helium temperatures due to the small value of T>), we find F a oc cos2<£> leading 
to 

vn? n £>A 2 L $lhv F /T~ 



12 eT 



12 A 



sin2y?. (12) 



A pure second harmonic is seen to develop, in complete agreement with the symmetry 
analysis. Note that in this temperature range the Meissner contribution to I ((f) becomes 
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Figure 6: (a) Meissner energy F M of an asymmetric 45° junction (in arbitrary units) 
as a function of the phase difference <p across the junction, (b) Superconducting current 
/ = (2n / & )dF M / dip (in arbitrary units) as a function of <p. In both figures a = 1, 2, 5, 10 
from top to bottom. 



negligible. This is a result of the thermal fluctuations which reduce the surface current 
Jo severely. 

Symmetric 45° junctions, i.e. junctions with 6± = —02 = 22.5°. In this case, trajec- 
tories with impact angles in the hatched regions of Fig. ^ are of type (i), whereas the 
remaining ones are of type (iii). As a function of the impact angle, the gap functions can 
be written Ai(tf) = A 2 (i?) = Acos(2$ - vr/4) and A 3 (V) = A 4 (tf) = Acos(2$ + tt/4). 
Again we have to make an assumption about the dependence of the transparency on the 
impact angle. We start by considering the case of a weakly dependent D(i9). As a rough 
approximation, we approximate all Andreev levels of type (iii) by their value for d = 0° 



2~Ifl ) and all levels of type (i) by the $ = 45° case, 



e(0,y>) = ±2- 1/2 A v /l-/J(0)sin 2 (<^/2) 



e(ir/4,(p) = ±2- 1/2 Asm( V /2)^/D(7i/4:), 

where we assume that D(tt/A) < D(0) 1. Note that for every impact angle $ there 
exist two levels with energies ±e(i?). Therefore J2 n = 1 an d the integral Eq. || 
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vanishes. In other words, the total supercurrent along the interface vanishes by symmetry. 
This does not mean, however, that no currents are flowing along the interface. In order 
to estimate the Meissner contribution to the interface free energy F{f, we model the 
current distribution generated by the Andreev levels as j(x) = Jq[5(x + £) — 8(x — £)] 
where £ is the coherence length and Jo = dxj y (x). The resulting magnetic field and 
the Meissner screening currents can be calculated within the London theory in a similar 
way as for the asymmetric junctions (see Fig. |5]) and the resulting contribution to the 
junction free energy is F^f = liqJq^. Note that F^f is reduced with respect to the case 
of asymmetric junctions (if we assume the same value of Jo in both cases) by a factor 
4£/A <C 1, since it arises from a higher order moment of the current distribution j y (x). 
Due to this small factor, F^ 1 /F A ~ (£/A) 2 //J even if we take the maximal possible value 
of Jo, Jo ~ ekpVp/d. Therefore the Meissner contribution to the interface free energy will 
be neglected and the Josephson current is calculated from 



27T 8F A kr? W 2 Be 
m = = i^E / ^oos^WJp. (13) 



This implies that j((p) = jn(y?) + ja{<p), where the normal and anomalous contributions 
are, respectively: 



n k F A /J(0)siny? Ay/1 - D(0) sin%/2) 

3n[if) ~ 8W 2 3/ J Jl-i)(0)8hi J W ' 



3a{ V ) = -f^ 008(^/2)^(^/4) tanh >- . (14) 



At intermediate temperatures A J D(tt/A) <C T <C A the anomalous (normal) Andreev 
levels are in the high (low) temperature limit. Therefore Eq. [14] simplifies considerably 
and the first two harmonics of j((p) can be written as 

- = D{0)-£d(*/4), 
Jl 2T 

l DW-h(^] D^/Af, (15) 



j L 8 w 24 \2Tj 

where jl = /l§y/2)kFA/Q Q d is comparable with the bulk depairing current den- 
sity. Thus, for weakly angle-dependent tunneling and above the energy of anomalous 
Andreev levels, theory predicts a sign change of the first harmonic at a temperature 
T* = AD(7r/4)/2D(0) and a second harmonic monotonically increasing with decreasing 
temperature. 

In the opposite limit of a peaked barrier transparency D($) = D(0) exp(— d/do), the 
junction is equivalent to a symmetric junction between s-wave superconductors with order 
parameters A/ a/2. Therefore the current-phase relation in the T — > and T — > T c limits, 
respectively, reads 



7T A 

jUp)R n = — 1= — 



D(0) 

sin ip sin 2<p 

16 
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Y 48 \2TJ y 



(16) 



Unlike Eq. 15, the result Eq. O predicts a monotonically increasing first harmonic with 



decreasing temperature. 

2.2. c-AXIS JUNCTIONS. Let us briefly discuss the supercurrent in c-axis junctions 
between YE^C^O^ (YBCO) and low-T c superconductors ||. In YBCO the dominant 
component of the superconducting order parameter has rf-wave symmetry 0. However, 
due to the orthorhombic structure of YBCO, a finite component with s-wave symmetry 
is admixed to the dominant rf-wave order parameter. The in-plane phase sensitive ex- 
periments imply that the (i-wave component remains coherent through the whole sample 
0, while an elegant c-axis tunneling experiment shows directly that the s-wave order 
parameter component does change sign across the twin boundary |26| . 

The above picture of the YBCO pairing state is challenged by the experimental ob- 
servation of a finite c-axis Josephson current between heavily twinned YBCO and a Pb 
counterelectrode ^7|. Namely, the contribution of the s-wave part of the YBCO order 
parameter to the Josephson coupling between a conventional superconductor (supercon- 
ductor with a pure s-wave symmetry of the order parameter, for example Pb or Nb) and 
YBCO should average to zero for equal abundances of the two types of twins in YBCO. 
In other words, the macroscopic pairing symmetry of twinned YBCO samples should be 
a pure (i-wave Q. Tanaka has shown that a finite second order Josephson current ob- 
tains for a junction between the s-wave and c-axis oriented pure (i-wave superconductors 



nfl . However, measurements of microwave induced steps at multiples of hf/2e on the 



I-V curves of Pb/Ag/YBCO tunnel junctions imply dominant first order tunneling [ 28 
Therefore the finite c-axis Josephson current has to result from a nonvanishing admixture 
of the s-wave component to the macroscopic order parameter of YBCO . Two alter- 
natives of how this can take place in the junctions based on twinned YBCO have been 
discussed in the literature: 

(i) Sigrist et al. have suggested that the phase of the s-wave component in YBCO 
does not simply jump from to ir upon crossing the twin boundary, but rather changes 



in a smooth way, attaining the value of tc/2 right at the twin boundary [29[. The twinned 
YBCO sample is thus assumed to exhibit a macroscopic d + is pairing symmetry. A 
related picture has been proposed by Haslinger and Joynt, who suggest a d + is surface 



state of YBCO W\ 



(ii) A difference in the abundances of the two types of twins implies a d + s symmetry 
of the macroscopic pairing state pi|. It has been pointed out p[ that also structural 



peculiarities of other type (such as a lamellar structure in a preferred direction) may lead 
to the d + s macroscopic pairing symmetry. 

The primary motivation of our experiment || was to determine which of the above 
two scenaria is realized in YBCO. In what follows we will show that this question can be 
conveniently answered by measuring the relative sign of the first two harmonics of I((p)- 
Let us consider first the d+ s scenario. As will be shown explicitly in Eqs. [H], 0, in this 
case the junction free energy can be written 

F(ip) = - F l cos ip + F 2 cos 2<p (17) 
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with Fx, F2 > 0. As shown in Fig. [7], this choice of signs corresponds to a free energy which 
is flat (curved) close to the minima (maxima). In terms of the current phase relation, 
this case corresponds to the usual case when !(</?) is steep in the vicinity of ±7r, as in the 
Kulik-Omelyanchuk theory [II]. 
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Figure 7: (a,b): Schematic drawing of the contributions from the first (solid lines) and 
second (dashed lines) harmonics to the free energy of c-axis junctions between a conven- 
tional superconductor and YBCO for various macroscopic pairing states in YBCO. (a): 
d + s pairing, (b): d + is pairing. (c,d): Total free energy (c) and the current-phase 
relation (d) in the d + s case (solid lines) and d + is case (dashed lines). 

The case of d + is pairing is more complicated. In fact, first order coupling between 
the conventional superconductor and the s-wave component leads to a term in free energy 
—Fx cos(v?conv — fs), where (p conv and <p s are the phases of the conventional superconductor 
and of the s-wave component in YBCO and F\ > 0. On the other hand, the second order 
term is the same as in the d + s scenario. Therefore, since <p = ^ conv — % where (fid is the 
phase of the d-wave component in YBCO, and since ip s — ipd = vr/2 in the d + is scenario, 
the junction free energy can be written 

F(ip) = -F 1 cos(tp - tp min ) - F 2 cos2(v9 - v9 min ), (18) 

where (p m i n = n/ 2 is the equilibrium phase difference and Fi 1 F 2 > 0, see Fig. ^. In the 
experiment it is a delicate issue to determine the absolute value of <p m i n . Usually f m i n is 
set equal to zero in the free energy minimum and with such conventions the two cases 
Eqs. p7T|JT~8| lead to a different relative sign of I\ and I2 for the d + s and d + is scenaria, 
respectively. 
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2.3. FACETED JUNCTIONS. In situations when the first harmonic l\ vanishes by sym- 
metry (i.e. for asymmetric 45° grain boundary junctions or for c-axis tunneling between 
an s-wave and a (i-wave superconductor), small local deviations from ideal geometry lead 
to a finite local coupling bl\ across the junction. In such cases the junction can be viewed 
as a parallel series of Josephson junctions with local critical current whose magnitude and 
sign fluctuate along the interface. It is customary to call the junctions with a positive 
(negative) bl\ as and tc junctions, respectively. Let us note in passing that very recently 
it has become possible to prepare samples with a prescribed pattern of and tc junctions, 
making use of the so-called zigzag junctions [j32fl . 

An example of a series of and 7r junctions is shown in Fig. where the geometry 
of an asymmetric 45° grain boundary between two twinned YBCO samples is shown. 
Twinning in the grain with Q\ = 0° is not shown explicitly, since it is irrelevant for our 
argument. Another often considered source of deviations from ideal interface geometry is 
associated with the experimentally well documented meandering of interface [|l^ . 




Figure 8: Schematic drawing of an asymmetric 45° grain boundary junction between 
twinned YBCO superconducting grains. In the upper grain the sign of the dominant lobe 
of the order parameter changes between the neighboring twins, leading to a sequence of 
and 7r junctions. Also shown are the spontaneously generated currents along the interface. 



If the average interface geometry is close to ideal, the number of and it junctions 
is roughly the same. As noted first by Millis |14| (for an alternative formulation making 
use of the sine-Gordon equation, see e.g. [T2|, [£|), this does not imply a vanishing 
Josephson coupling across the interface. The general idea is as follows: let the average 
phase difference across the junction is <p. The phase difference along the interface is 
modulated from this value towards and ir in the and it junction regions, respectively, 
by an amount ~ \. This leads to spontaneously generated currents along the interface 
which are schematically depicted in Fig. ||. In this way, the junction gains Josephson 
energy cx \i whereas a Meissner energy of only oc x 2 is lost. The total energy gain 
is maximal for ip = ±7r/2 and this explains qualitatively the development of a second 
harmonic in such junctions. 

In what follows, let us discuss the faceted scenario in more detail. For definiteness, we 
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Figure 9: Schematic drawing of a c-axis junction between an s-wave superconductor and 
a twinned YBCO sample. The black and white regions (of typical size a) correspond to 
and it junctions, respectively. Also shown are the penetration depths in the s-wave 
superconductor (X R ) and in YBCO (A, A c ). 



consider the case of a c-axis junction between an s-wave superconductor and a twinned 
YBCO sample, see Fig. [| For the sake of simplicity, we model the randomly twinned 
structure by a checkerboard distribution of and ir junctions with a periodic distribution 
of the local critical current density j c ( x ) — joQ^), where g(x) = sin(7rx/a) sin(7n//a) and 
for a we take the typical size of the twins. Here we have assumed a coordinate system 
such that the interface lies in the xy plane. According to |14]], a spontaneous magnetic 
field is generated, which we now find by minimization of the total junction energy as a 
function of the phase difference <p across the junction. The symmetry of the problem 
dictates that the magnetic field reads 



B x = B n exp(—K n \z\) sm(mix/a) cos(mry / a) , 



By = — B n exp(—K n \z\) cos(mix/a) sm(mry/c 



B, 



0. 



(19) 



where the prime restricts the summation to odd values of n. Note that V • B = 0, as 
it should be. For z > 0, i.e. in the conventional superconductor, the B field satisfies 
the equation (V 2 - A^ 2 )B = if K 2 n 

the anisotropic cuprates, B 



X" 2 



1 + 2n 2 (tcXr/clY 
(B x , B y , 0) is governed by B x>y = 



On the other hand, in 



B XtV 



1 + 2n 2 (vrA c /a) 2 . From Eq. |T|, the 
culated using VxB = /x j. The total 



and therefore for z < we have to take K 2 = A -2 
distribution of Meissner screening currents can be ca 

junction energy per unit area corresponding to the magnetic field described by Eq. [19] is 

dSg(x) cos 6(x), (20) 



4/i n 

where S is the junction and 



1 Vr 2 a $oJo 



2nS 



X n = X R ^/l + 2n\Tx\ R /af + Ayi + 2n 2 (rr\ c /a) 
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The local phase difference across the junction 8(x) can be determined from the Josephson 
equations 



86 

dx 

m 

dy 



^ \X 2 j x (z = 0_) - XIjAz = + )l = ^ £ B n A r 

27T/i 



A 2 j y ( 2; = 0_)-A^(z = + ) 



cos sm , 

a a 

niry 



2n J-^ . nirx 

— > B n \ n sm cos 

$ n a a 



(21) 



These equations are solved by 6>(x) 
and 



ip + x( x ); where (p is the average phase difference 



, . 2a ^ -B n A n . nnx . niry 

y(x) = — > sm sm 

U $ „ « a a 



(22) 



is the spontaneously generated modulation of #(x). 

The Fourier components B n have to be determined by minimization of the energy 
Eq. To this end let us assume now that \x\ "C 1, which assumption will be justi- 
fied at the end of the calculation. The second term in Eq. ^ simplifies in this case to 
(2/i ) _1 i?i-BcfrAi sin ip, where B e{i = (2n)^ 1 fi j a. Note that Fourier components B n with 
n > 3 raise the energy and therefore vanish. This means that the junction energy is at the 
end a function only of B\, E n = (Ai/4/x ) (-£>i + B eS sin tp) 2 — B 2 S sin 2 ip . The minimum 
is reached for B\ = —B c g simp. The minimized value of the junction energy density and 
the resulting current-phase relation read 



2 / u $o 



sin 2 ip, 



sin 2<p. 



(23) 



Note that j oc j 2 oc T> 2 . Therefore the second harmonic scales with the barrier trans- 
parency in the same way both within the flat (for T ^> T>A) and the faceted scenaria. 
Finally, let us point out that the result Eq. ^3| is valid only for \\\ <C 1, or equivalently 
for 

(24) 



2.4. THE JOSEPHSON PRODUCT. According to standard theory (for homogeneous 
featureless barriers), the Josephson product I\Rn (where Rn is the junction resistance 
in the normal state) is independent of the junction area and of the barrier transparency, 
thus giving an intrinsic information about the superconducting banks. The measured 
Josephson product of cuprate grain boundary junctions |7] can be well described by 
IiRn = a 2 (7r/4)(A/e) COS26 1 , where A is the maximal superconducting gap. The angular 
dependence of IiRn is fully consistent (apart from the ^-independent renormalization 
factor a 2 ~ 1CT 1 ) with the BCS prediction Eq. |] for rough junctions between <i-wave 
superconductors. This means that 

(I c + I s )R N = a 2 (n/2)(A/e), (25) 
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much smaller than the theoretical prediction Eq. [3] for I c Rn- Note that the interpretation 
of a small factor a 2 as a result of a nearly complete cancellation of I c and I s is not 
plausible, since the cancellation would have to occur for all misorientation angles, whereas 
the physical origin of the I c and I s terms is quite different. 

It is worth pointing out that Eq. ^ applies (with the same a 2 ~ lO^ 1 ) also for the 
break junctions (in which the misorientation angle 8 = 0) p?3| . Moreover, in [3~3| it has 



been shown that A is not depressed in the junction region, thus explicitly demonstrating 
the breakdown of the BCS prediction for I\Rn in the cuprates. 

There exists no generally accepted explanation of the small renormalization factor 
a 2 . One of the reasons is that the microstructure of Josephson junctions is typically 
quite complicated. In fact, it is well known that small angle grain boundaries can be 
modelled by a sequence of edge dislocations, while at larger misorientation angles the 
dislocation cores start to overlap and no universal picture applies to the structure of 
the grain boundary. For large-angle grain boundaries, Halbritter has proposed []35| that 
the junction can be thought of as a nearly impenetrable barrier with randomly placed 
highly conductive channels across it. If due to strong Coulomb repulsion only the normal 
current (and no supercurrent) is supported by these channels, the small value of I\Rn 
follows quite naturally. 

On the other hand, we have argued |36|] that the smallness of the Josephson product 
does not follow from the particular properties of the barrier, but is rather an intrinsic 
property of the cuprates. Such a point of view has been first advocated in |]37| . However, 



that paper did not consider alternative more conventional explanations. In order to 
support our point of view, in [^] we have discussed the Josephson product for intrinsic 
Josephson junctions in the c-axis direction. Such junctions can be viewed as an analogue 
of a6-plane break junctions (since the misorientation angle vanishes for both), but are 
preferable because of simpler geometry of the interface. Moreover, zero energy surface 
bound states which may develop at afe-plane surfaces because of the (i-wave symmetry 
of the pairing state [0] do not form in the c-axis direction, simplifying the analysis of 



intrinsic Josephson junctions. By an analysis of the experimental data [38|, we came to the 



conclusion that the experimental Josephson product is reduced with respect to theoretical 
predictions. Thus we have shown that although the barriers in grain boundaries and in 
intrinsic Josephson junctions are of very different nature, both types of junctions exhibit 
a suppressed Josephson product. Therefore we concluded that this suppression is not 



due to specific barrier properties as suggested in pSj , but rather due to some intrinsic 
property of the high-T c superconductors. 



3. EXPERIMENT 

As shown in the previous Section, theory predicts anomalous behavior of the Josephson 
current in the cuprates as a function of the phase difference (p and temperature. One of 
the most striking predicted phenomena is a non-monotonic temperature dependence of 
the Josephson critical current, forced by the negative contribution of anomalous Andreev 
states. However, this has not been observed by standard means in spite of enormous efforts 
of many groups. The main obstacle seems to originate from the interface roughness which 



17 



leads to a suppression of the anomalous Andreev states [10] . Even the most promising 
type of Josephson junctions, grain boundary junctions, contain defects on a characteristic 



scale of 1 /iin (Fig. 2 of [12|) coming from defects in the bicrystal substrate, see Fig. |10 



A similar problem is relevant also in c-axis Josephson junctions, where a large second 
harmonic of the Josephson current is expected. In that case, in order to avoid tunneling 
in the ab plane direction and to achieve a high transparency of the interface, atomically 
flat surfaces are required. 



bicrystal boundary- 




Figure 10: Schematic picture of the RF SQUID (taken from |4|]). The YBCO thin film 
occupies the gray area. The inset shows an electron microscope image of the narrow grain 
boundary Josephson junction. 

Thus, in order to observe the predicted anomalies, one should avoid the structural 
defects at the interface. Following this strategy the c-axis Josephson junctions were pre- 
pared in situ by covering YBCO by thick Au layer preventing the degradation of the 
YBCO surface during processing. On the other hand, we have chosen the position of the 
grain boundary junctions in such a way that they were placed between the defects of the 
substrate shown in Fig. ID. Therefore we were forced to reduce the size of the junction 
below 1 /im. However, for submicron 45° junctions where the anomalies should be most 
pronounced, the Josephson coupling energy Ej is comparable with T in the temperature 
range of interest. This leads to large phase fluctuations and consequently to a suppression 
of the Josephson current. At first sight it seems that the Josephson current of submicron 
45° junctions cannot be measured. But this is not so. 

3.1. EXPERIMENTAL METHOD. Phase fluctuations can be suppressed by connecting 
the Josephson junction into a ring (rf SQUID), since in that geometry the phase differ- 
ence across the junction is controlled by the large phase stiffness of the superconducting 
ring. This technique therefore offers a unique possibility to study Josephson junctions at 
temperatures T much higher than the junction energy. In addition, one can change the 
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phase difference f on the Josephson junction by applying external magnetic flux to 
the ring 

2ttL s 

< P = ( Pdc-—f—I{<P), (26) 

where fdc = 27r&dc/&o, L s is the inductance of the rf SQUID, and 1(f) is the Josephson 
current. We have used the modified Rifkin-Deaver method |4l], f|2] to restore the current- 
phase dependence 1(f). The method is simple: The rf SQUID is inductively coupled to 
a high-quality parallel resonance circuit driven at its resonant frequency ujq. The angular 
phase shift a between the rf driving current J r / and the voltage across the circuit is 
measured by a rf lock-in voltmeter as a function of the external magnetic flux. The 
flux is set by a dc current I<f c through the coil of the resonant circuit with inductance 
Lt- Thus the total external magnetic flux can be expressed as f e = fd c + frf, where 
fdc = ^irMIdc/^o and f r f = 2irMI r f /$o and M is the mutual inductance between the rf 
SQUID and the resonant circuit. The rf SQUID and the resonant circuit are characterized 
by quality factors q = u L s /R and Q = Rt/uj L t , respectively. The analysis of such a 
system can be considerably simplified in the adiabatic (g < 1, Q > 1), small signal 
regime (f r f C 1 + (2TcL s /^ )dI(f)/df), when the rf SQUID follows adiabatically the 
signal from the resonant circuit and the internal flux can be expressed by Eq. [26|. I(fd c ) 
is calculated from the experimental a(fd c ) data using 

I i ( fdc) = 7^-;- / tan a(f dc )d f dc , (27) 
27tQ<P Jo 

where Jo is the period of a(Idc) (i-e. MI = $ ) and the quality factor Q is measured 
independently from the width of the resonance curve of the parallel resonant circuit. 
Using Eqs. |26|j27| 1(f) can be restored. This method, being differential with respect to f, 
provides a high sensitivity of the current phase measurement. 

Let us emphasize that the critical current determined by this method does not depend 
on the inductance of the rf SQUID. In fact, one can easily show that I c = I(f dc ), where 
tan a(f dc ) — can be determined from experimental data by requiring J n tan a(fd c )d fdc = 
0. The last condition comes from the periodicity of 1(f) and enables a subtraction of a 
constant phase shift coming from electronics and cables. Since the only sample charac- 
teristics entering Eq. [27|is 1(f), the measurement system can be calibrated using samples 
with known critical current and 1(f)- We have usually carried out the calibration making 
use of Nb rf SQUIDs. Typical values of the quantities used in our experimental setup are: 
L s = 80 pH, L t = 0.27 /iH, Q = 155. Critical currents measured by the present method 
and by standard transport measurements were in good coincidence. As already explained, 
1(f) is measurable even if the thermal energy exceeds the Josephson coupling energy. In 
fact, critical currents down to 50 nA were recently detected at T = 4.2 K [j|] using cooled 
preamplifier placed near the parallel resonant circuit. All 1(f) measurements presented in 
this paper have been performed in a gas-flow cryostat with a five-layer magnetic shielding 
in the temperature range 1.6 < T < 90 K. 

3.2. GRAIN BOUNDARY JOSEPHSON JUNCTIONS. The rf SQUIDs were prepared 
at IPHT Jena by laser deposition of YBCO thin films of thickness 100 nm on bicrystal 
substrates. They were patterned in the shape of a square washer 3500 /mix 3500 /im with a 
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hole 50 /mix 50 /iin by electron beam lithography (see Fig. ITU). The estimated Josephson 
penetration depth Aj is much smaller than the width of the wide junction, wi = 1725 
fim, and larger than the width of the narrow junction, w s = 1 /im and w s = 0.7 /im for 
asymmetric and symmetric 45° junctions, respectively. Thus the behavior of the rf SQUID 
is dictated by the narrow junction only. For symmetric 45° junctions the submicron bridge 
was formed at a position between the defects of the substrate which are visible in Fig. |T0[ 
The experimental results are shown in Fig. [11]. Local minima appear at low temperatures 
on the a((pd c ) curve close to (fd c = 27m where n is integer. The existence of the local 
minima dictates 



d(p 3 



> 0. 



(28) 



ip=2irn 



Note that neither the conventional tunneling theory nor the I and II theories by Kulik 
and Omelyanchuk predict such local minima on the derivatives of CPR M. 
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Figure 11: The phase angle a as a function of (pd c measured at different temperatures, 
(a) Asymmetric 45° grain boundary ||. Top to bottom curves correspond to T =40, 30, 
20, 15, 10, and 4.2 K, respectively, (b) Symmetric 45° grain boundary [|J. Top to bottom 
curves correspond to T =35, 30, 25, 20, 15, 11, 10, 5, and 1.6 K, respectively. 

To demonstrate that the anomalous behavior of the Josephson current is a peculiarity 
of 45° junctions we have measured /(<£>) of symmetric YBCO grain boundary junctions 
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with misorientations 9 = 24° and 36°. The results of this study are shown in Fig. 0. We 
have found that the critical current decreases exponentially with the misorientation angle 
9, exp(—9/9o) with 9 = 5.6°, in good agreement with [)43[| . Moreover, both for 9 = 24° and 
36°, no considerable deviation from sinusoidal dependence of I(ip) was observed, implying 
tunneling as the dominant transport mechanism. 




1.00 



Figure 12: Current-phase relation for asymmetric (a) and symmetric (b) 45° grain bound- 
ary junctions (taken from H). (a) Top to bottom curves at <p = 0.757T: T = 4.2, 10, 15, 
20, and 30 K. (b) Bottom to top curves at <p = 0.375vr: T = 1.6, 5, 10, 11, 15, and 25 K. 

Let us return to the case of 45° junctions now. Theory predicts that the first harmonic 
should be suppressed for such junctions. If we in addition take into account that the 
barrier transmission is low, we can model 1(f) of such junctions by neglecting terms 
of higher orders (n > 2) in Eq. |l|. Therefore the condition (^) for the existence of 
the local minima at (fdc = 27m dictates < —1/8. Thus we conclude that the 45° 

junctions exhibit anomalously large second harmonic of the Josephson current which has 
opposite sign in comparison with I\. The current-phase relations calculated from the 
experimental data in Fig. [11] are shown in Fig. |T^. Note the anomalous form of I ((f) 
at low temperatures. The temperature dependence of the first two harmonics I\ and I2, 
determined from a Fourier analysis of I [(f), is shown in Fig. [13]. For both asymmetric and 
symmetric 45° junctions the second harmonic monotonically increases as the temperature 
decreases. The first harmonic is more or less constant for asymmetric junctions whereas for 
symmetric 45° junctions the first harmonic starts to exhibit a downturn below 20 K. The 
most striking result is that for T=12 K, Ii changes sign. In the same temperature region 
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where I\ starts to exhibit a downturn, the absolute value of I2 rises from a negligible value 
at high temperatures to values comparable to I\ at low temperature. This experimental 
fact suggests a common origin of both phenomena. Similar effects were theoretically 
predicted for Josephson junctions between <i-wave superconductors (Section 2) and our 
experiments can be therefore regarded as an independent experimental test of the <i-wave 
symmetry of pairing in YBCO. 
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Figure 13: Critical current I c (triangles) and the first two harmonics I\ (squares) and 
I2 (circles) as a function of temperature for asymmetric (a) and symmetric (b) grain 
boundary 45° junctions (taken from ||). The lines are guides to the eye. 



3.2. c-AXIS JOSEPHSON JUNCTIONS. The c-axis Josephson junctions were fabricated 



and characterized jointly by the IRE Moscow and Chalmers groups following [44]. Epi- 
taxial (OOl)-oriented YBCO thin films with thickness 150 nm were obtained by laser de- 
position on (100) LaAlOs and (100) SrTiOs substrates and in situ covered by a 8^-20 nm 
thick Au layer, thus preventing the degradation of the YBCO surface during processing. 
Afterwards, 200 nm thick Nb counterelectrodes were deposited by DC-magnetron sput- 
tering. Junctions with dimensions 10 x 10 fim 2 were formed by photolithography and low 
energy ion milling techniques. The interface resistance per unit area Rn = RnS (where 
Rn is the normal state resistance and S is the junction area) was Rn = 10 5 ^- 10 6 Q-cm 2 . 
Details of the junction fabrication were reported elsewhere |44]| . 

Surface quality of the YBCO films is very important when current transport in the c- 
direction is investigated. High-resolution atomic force microscopy reveals a smooth surface 
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consisting of approximately 100 nm long islands with vertical peak-to-valley distance of 
3-h4 nm ||. We can exclude that substantial afr-plane tunnel currents flow between YBCO 
and Nb at the boundaries of these islands. In fact, theory predicts formation of midgap 
states at the surface of semi- infinite Cu0 2 planes [|13|, £EJ. Therefore zero bias conductance 
peaks should be expected in the I-V characteristics at temperatures larger than the critical 
temperature of Nb, if the contribution of a6-plane tunneling was nonnegligible. However, 
no such peaks have been observed for all fabricated Nb/Au/YBCO junctions. Moreover, 
from the size of the islands and from the vertical peak-to-valley distance we estimate 
that the area across which afr-plane tunneling might take place is only 6% of the total 
junction area. However, since the interface resistances per square are of the same order 



of magnitude |27| for both, the c-axis and afr-plane junctions, we conclude that afe-plane 
tunneling from YBCO, if present, is negligibly small. 

More than 20 junctions were characterized by transport measurements. At small 
voltages typical I-V curves can be described by the resistively shunted junction model 
with a small capacitance Typical critical current densities were j c — l-j-12 A/cm 2 and 
j c Ra = 10 -T- 90 //V. The differential resistance vs. voltage dependence Rd{V) exhibits a 
gap-like structure at V ~ 1.2 mV. This structure has a BCS-like temperature dependence 
and disappears at T c r « 9.1 K, therefore we ascribe it to the superconducting energy gap 
of Nb. 
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Figure 14: Phase shift a as a function of ipdc for a c-axis junction at T—1.7, 2.5, 3.5, 4.2, 
and 6.0 K (from bottom to top). Taken from ||. 



The current phase relation of the c-axis Josephson junctions was measured by clos- 
ing the Nb/Au/(001)YBCO heterostructure into a superconducting ring with the same 
geometry as for 45° grain boundary junctions. The experimental results are shown in 
Figs. When compared with 45° junctions, the second harmonic of the Josephson 

current was considerably smaller but still anomalously large, leading to local minima of 
a((fdc)- As follows from the analysis in Section 2.2, the opposite signs of the first and 
second harmonics of the Josephson current provide direct evidence that in our YBCO 
samples, pairing with a macroscopic d + s symmetry is realized. The large first harmonic 
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has to be due to an uncompensated s-wave component whose origin is discussed in detail 
in the next Section. 
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Figure 15: The current-phase relation I(<p) of the c-axis junction from Fig. [14] at 
T=1.7, 2.5, 3.5, 4.2, and 6.0 K (from top to bottom). Inset: Temperature depen- 
dence of I\ (squares) and (circles). Solid lines are fits to Eqs. ([H],^) using 
Afl(T) = A fl (0) ta,nh[A R {T)T cR /A R {0)T]. Taken from [|. 



4. DISCUSSION 

4.1. GRAIN BOUNDARY JUNCTIONS. Asymmetric 45° junctions. The large second 
harmonic observed in || confirms the naive expectations based on the symmetry of the 
pairing state in the cuprates. However, at the time of writing the paper ||, we were not 
able to explain the details of the shape of the current-phase relation. In particular, we 
could not find any mechanism leading to a current-phase relation which was steep in the 
minima and flat in the maxima of energy. For instance the Kulik-Omelyanchuk theory 
leads to an exactly opposite picture of jumps in I(<p), which occur in the maxima of 
energy (at ip = ±7r). This is a consequence of the fact that for ip = ±7r two locally stable 
branches of the junction energy cross, one centered around <p = and another one around 
ip = ±2tt. We believe that in the present work we have found one physically plausible 
scenario for how the pattern of jumps of I(ip) observed in || can be explained. In fact, 
Fig. |6| shows that behavior qualitatively similar to the results of || obtains if the Meissner 



energy dominates the junction energy Eq. 11. A more detailed investigation of this idea 



is under way. A preliminary analysis indicates that in order to be applicable to [[|, also 
the effect of a spatially fluctuating barrier similar to that discussed at the end of Section 
4.2 needs to be taken into account. 

Although consistent with naive expectations, under closer inspection, the experimen- 
tally found large value of the second harmonic remains mysterious: 
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In the flat scenario, the second harmonic is given by Eq. |T2|. On the other hand, 
by combining Eqs. [4|,p5| we estimate that the first harmonic jiRa ~ ira 2 5A/e, where S 
is the deviation of the misorientation angle 9 from the nominal value of 45°. Taking a 
realistic value of such a deviation at least 5 ~ 10 -2 , by comparing to Eq. [12] we find that 
the experimental value J2/J1 ~ 1 can be explained, only if i? ~ 1 (and therefore the 
tunneling has to be possible in a wide range of impact angles). In this estimate we have 
made use of the well established []7|] value T> ~ 10~ 4 and we have taken A/T ~ 10 2 at 
helium temperatures. Note that in order to obtain J2/J1 ~ 1 we had to assume that only 
the first harmonic is renormalized by a 2 ~ 1CT 1 , while the second harmonic has to remain 
unrenormalized. This implies very peculiar microscopic physics. 

Now let us consider the faceted scenario. As an order of magnitude estimate, we use 
the result Eq. |23| for the second harmonic. Unlike in the c-axis case, in the present case all 
penetration depths have to be set equal to the in-plane penetration depth of the cuprates 
A ~ 0.15yum. This yields j 2 = (Sk&q)' 1 noj 2 a 2 \i, where Ai = 2Ay / l + 2(7rA/a) 2 . For a we 
take the typical length scale of faceting, a w 0.1/im, since the twin size is even smaller. 
This yields Ai ~ 2/im. The only unknown in the expression for j 2 is the local current 
density j - 111 what follows we determine j from j = y / 87rj2 < £o/(A i o a2 ^i) where for j 2 
we take the experimental critical current density j 2 ~ 10 4 A/cm 2 0, and we obtain a 
reasonable value jo ~ 10 6 A/cm 2 . Note that since _B e fr = (27r) _1 /iojo a ~ 3 x 10~ 4 T 
and $o/(2«Ai) ~ 5 x 10~ 3 T, the criterion Eq. |24] for the applicability of Eq. [23] is well 
satified. On the other hand, we estimate the first harmonic from ji/jo ~ (o-fL) 1 ^ 2 , which 
is a random walk-type formula, indicating that ji averages to zero in a sufficiently long 
junction. The experiment requires j 2 > ji- This is possible only for sufficiently long 
junctions, L ~ 1 mm, whereas in || much shorter junctions with L ~ 1/zm were studied. 
Therefore within the faceted scenario it is not possible to explain the large measured I2/I1 
ratio. 

We conclude that both within the flat and faceted scenaria, the second harmonic 
measured in appears to be anomalously large. 

Symmetric 45° junctions. The anomalous temperature dependence of the first har- 
monic and the large second harmonic observed in this type of junctions is qualitatively 
consistent with the theoretical prediction Eq. [15| and not with Eq. IB]. The data seems to 
imply (as was the case for asymmetric junctions as well) that the impact angle dependence 
of the barrier transmission is weak. 

In what follows we attempt a more quantitative discussion of the results found in |4[] . 
According to Fig. 13 the first harmonic changes sign at a temperature T* w 12 K. This 
together with the estimate A « 20 meV [0 requires D(7r/4:)/D(0) ~ 0.1, i.e. the barrier, 
although presumably quite thin, can't be modelled by a delta function. 

Now let us compare the relative magnitudes of the first and second harmonics. Theory 
predicts that for T > T*, the maximal value of j\ is j^^/ji ~ D(0). On the other hand, 
from Eq. |T5]it follows that at T = T* the second harmonic \ j 2 {T*)\/j L ~ j D(0) 3 /24D(tt/4). 
Therefore according to theory |j2(r*)|/jf iax fa D (0) 2 / '24D (it / '4) ps D(0)/2, whereas from 
Fig. [H3]we estimate \j2{T*)\/jf l3 ' x ~ 1. Thus theory can describe the experimental results 
only if D(0) ~ 1. 

However, in what follows we show that D(0) C 1 and therefore the experimental sec- 
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ond harmonic is again too large when compared with simple minded theory. In fact, in or- 
der that higher-order harmonics are negligible |4j] even at the lowest studied temperatures 
T m i n ~ 1.6 K, we require that the typical mid-gap state energy Ay / D(7r/4)/2 < 2T min , 
yielding D(n/4) < 4 x 1(T 4 , and therefore D(0) « 10D(7r/4) < 4 x 10~ 3 . Moreover, 
D(0) ~ 10~ 3 seems to be consistent with the exponential decrease of j c with the misori- 
entation angle 9 [Q. 

In || we have proposed that the first harmonic might by suppressed by interface 
roughness according to Eq. f|. This would require a very small roughness parameter 
x ~ 10 -3 and therefore the junction would have to be nearly completely rough. However, 
we have overlooked the fact that for such junctions it is hard to believe that the deviation 
5 from the misorientation angle 9 = 45° can be sufficiently small in order to keep the first 
term in Eq. ^ small. 

In summary, the second harmonic of 45° grain boundary junctions seems to be too 
large to be explicable by conventional theory. 



4.2. c-AXIS JUNCTIONS. Let us start by estimating the transparency of the barrier 
between YBCO and Nb from the normal-state resistance per unit area Rn- According to 
the band-structure calculations (for a review, see [|B|), the hole Fermi surface of YBCO 
is a slightly warped barrel with an approximately circular in-plane cross-section (to be 
called Fermi line) with radius kp. In what follows, we represent the electron wavevector 
k in cylindrical coordinates, k = (k,8,k z ). We estimate the uncertainty of the in-plane 
momentum as 5k « 2tt/1, where I is the characteristic size of the islands on the YBCO 
surface (fig. 1) . We evaluate R u making use of the Landauer formula and note that only 
tunneling from a shell around the Fermi line with width 5k is kinematically allowed. 
The barrier transparency D(9) depends on the details of the c-axis charge dynamics in 
YBCO, with maxima in those directions 6, in which the YBCO c-axis Fermi velocity 
w(8) is maximal. Since for 9 = 7r/4 and symmetry equivalent directions w(9) is minimal 
|[47|| ) we expect that there will be 8 maxima of D(9) on the YBCO Fermi line where 
D(9) ~ D, which are situated at 9 = 9q and symmetry equivalent directions. The 
modulation of the function D(9) depends on the thickness of the barrier between YBCO 



and Nb fl48" |. We consider two limiting distributions of the barrier transparency D(9) 
along the YBCO Fermi line: (a) a featureless D(9) ~ D and (b) a strongly peaked D(9), 
roughly corresponding to thin and thick barriers, respectively []48"| |. In the thick barrier 



limit the angular size of the maxima of D(9) can be estimated as 59 « 5k/ kp. With these 
assumptions we find 

Rn 1 = (29) 

where A measures the number of conduction channels and (. . .) denotes an average over the 
junction area. In the thin and thick barrier limits, we find A « kp5k/n and A « 25k 2 /tt, 
respectively. Taking I w 100 nm and kp w 0.6 A -1 ||4"9| , the measured R n = 6 x 10~ 5 
ficm 2 can be fitted with (D) thin w 1.7 x 10~ 5 and (D) thick w 8.3 x 10" 4 . 

Now we can turn to the discussion of 1(f). Since we have observed no midgap surface 
states in the Rd(V) curves, we can neglect the surface roughness, and the Josephson 
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current can be calculated from IKO 



77 ) = ^\^k TV £>Aj;A k siny 

where the sum over k, 6 is taken over the same regions with areas A as in Eq. (|29|) , Ar 



and Ak are the Nb and YBCO gaps, respectively, and Vt,i = y uj 2 + Af with i — R, k. 
Keeping only terms up to second order in the (small) junction transparency D, the first 
and second harmonics of the Josephson current densities for an untwinned YBCO sample 
read 

■ ( r r\ d ~ A S A R (T) 

jo{T)R n « — — - — , (31) 
J 2 (T)i? D - --———tub^^), (32) 

where AJJ = irAd[2 ln(3.56Arf/T c ij)] _1 and AjJ = A^| cos2#o| in the thin and thick barrier 
limits, respectively. In Eqs. ([H],[!^) we have assumed that the angular variation of the 
YBCO gap can be described as A(9) = A s + A d cos2#, where A d and A s are the d- 
wave and s-wave gaps. We have assumed that A^ is larger than both, A R and A s . The 
factor A s /Aj can be estimated from the measured joRn products for Josephson junctions 
between untwinned YBCO single crystals and Pb counterelectrodes. For such junctions 
jo(0)R n « 0.5 -=- 1.6 mV [0. Using the Pb gap A R = 1.4 meV in Eq. we obtain 



A s /A^ w 0.36 -=-1.1. 

From the relative sign of I\ and I2 we know that the finite first harmonic has to be due 
to the macroscopic d + s symmetry of our YBCO sample. The simplest way how this can 
be realized is to assume that the numbers of the two types of twins are unequal. In fact, 
detailed structural studies show that this is the case for sufficiently thin YBCO films even 
if they are grown on the cubic substrate SrTiOs [ETJ] . If we denote the twin fractions as 
(1 +8)/2 and (1 — S)/2, then the measured first harmonic of the CPR, ji, is proportional 
to the deviation from equal population of twins, j\ = 5j |TL|. Using A# = 1.2 meV 
determined from the Rd(V) data and our estimate A S /AJ ~ 0.36 -r 1.1, we find that the 
measured first harmonic j\ can be fitted with 5 ~ 0.07 -7- 0.21, which is in qualitative 



agreement with [51], where S ~ 0.14 for 1000 A thick YBCO films has been observed. 

Fitting the measured j2-Rn by Eq. (0), (D 2 ) / (D) m 3.2 x 10 -2 was found, which is 
much larger than both (-D)thin and (-D)thick- Similarly as in the case of grain boundary 
junctions, this means that the experimental second harmonic is too large when compared 
with naive estimates. In what follows we discuss the possible causes of such behavior. 

Let us start by considering the faceted scenario. From the known values of joRa W? \ 
and R n we estimate jo ~ 8 -r- 27 A/cm 2 . Taking for a typical twin size a ~ 10 nm, the 
effective magnetic field from Section 2.3 is estimated as B e ^ w (1.6 -j- 5.4) x 10 -10 T. 
Moreover, since \ R ~ 39 nm, A ~ 240 nm, and A c ~ 3 /im (as an upper bound, for 
the cuprates we take numbers which are valid for underdoped YBCO ||52||), we find Ai 



320 fim. Since this implies $o/(2aAi) ~ 3 x 10 4 T, the criterion Eq. [24] is seen to be 
well satisfied and therefore the second harmonic can be calculated from Eq. E3|, yielding 
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1.72 1 -Rn ~ (3 4- 35) x 10 11 V, orders of magnitude smaller than the experimental value 
\j 2 \Ru « 15 fiV. 

Thus we were led to look for alternative explanations of the large second harmonic. 
In HH we came up with an explanation which assumed that the junction transparency D 
is a fluctuating function of the position r. Adopting the WKB description of tunneling 
4~8"|, we assumed that the local barrier transparency D(s(r)) = exp(— s — s(r)), where Sq 



is the WKB tunneling exponent and s(r) its local deviation from the mean. Assuming a 
Gaussian distribution of s with a mean deviation rj, P(s) oc exp(— s 2 /r] 2 ), we estimated 
the spatial averages as (D n ) = dsP(s)D n (s). In the thin barrier limit, the values 
(£> 2 )thin = 8.6 x 10~ 7 and (-D)thin = 2.3 x 10~ 5 required to fit the experiments correspond 
to an average WKB exponent SQ hin ~ 15.5 with ?7 t hin ~ 4.3. In the thick barrier limit we 
obtain Sg hlck « 9.1 and ?7thick ~ 2.8. 

4.3. MICROSCOPIC IMPLICATIONS. Very recently it has been pointed out by one of us 
36fl that the two apparently unrelated experimental facts, namely the suppressed Joseph- 



son product I\Rn and the enhanced ratio \l2/h\, can be explained by a single assumption 
that in the cuprates some mechanism is operative which leads to a suppression of I\, while 
leaving and I2 intact. In what follows we describe one such mechanism which we be- 
lieve to be the most promising one. Namely, we suggest that at low temperatures the 
superconducting state of the cuprates supports fluctuations of the superconducting phase. 
Such fluctuations presumably do not affect R^, while they do influence the Josephson cur- 
rent. In simplest terms, if we denote the phases of the superconducting grains forming 
the junction as fa, then the fluctuations renormalize the first and second harmonics by 
the factors (e lr ^ 1 )(e 1 ^ 2 ) and (e 2 ^ 1 ) (e 2 *^ 2 ), respectively, where (. . .) denotes a ground-state 
expectation value. Thus experiment requires that the fluctuations have to be of such type 
that |(e 4<?i )| = a ~ 0.3 and |(e 24<?i )| ~ 1. Precisely this behavior is expected if the c?-wave 
order parameter fluctuates towards s-wave pairing (which pairing is expected to be locally 
stable within several microscopic models of the cuprates). 

An independent check of our picture is provided by measurements of the Joseph- 
son product for junctions between the cuprates and low-T c superconductors. In such 
a case, we predict that the Josephson product should be renormalized by a factor a 
instead of a 2 for junctions between two cuprates. Junctions of this type have been stud- 
ied extensively in the past. In particular, afe-plane junctions between YBCO and Pb 
exhibit Josephson products of 0.2 4- 1.2 mV p7| . On the other hand, Eq. |3] predicts 
I c Rn ~ e _1 A Pb ln(4A Y Bco/Apb) ~ 5.7 mV in this case, if we assume a (100) YBCO 
surface and take Aybco ~ 20 meV and Apb ~ 1.4 meV. We interpret the large experi- 
mental scatter of I c Rn as being due to varying afe-plane tunneling direction. Thus theory 
has to be compared with the largest experimental value, and the theoretical result has to 
be multiplied by a renormalization factor a ~ 0.2 in order to bring it in agreement with 
experiment. This value is in semiquantitative agreement with a ~ 0.3 which was deter- 
mined from the Josephson product of grain boundary junctions. It is worth mentioning 
that the phase fluctuation picture also may be relevant for the experiment ||, where a 
large second harmonic has been found in a c-axis Josephson junction between YBCO and 
Nb. 
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5. CONCLUSIONS 



In this paper we have tried to argue that the Josephson effect can provide nontrivial 
information not only about the symmetry of the pairing state in the cuprates, but also 
about the fluctuations of the superconducting order paramater. The latter are expected 
to be quite large especially in the underdoped region, where charge fluctuations should 
be suppressed. A preliminary analysis of the Josephson product and of the current-phase 
relation in grain boundary and c-axis junctions indicates J36| that the phase fluctuations 
of the superconducting order parameter are quite large and of a very special type, favoring 
local fluctuations from c?-wave towards extended-s pairing. 

Surprisingly, the potential of the Josephson effect as a test of the phase rigidity of 
the cuprates seems to have been missed in the past and therefore quite few studies have 
attempted quantitative analysis of the Josephson phenomena. Therefore there are still 
more questions than answers in this field. Some of the most pressing problems of the field 
(from our point of view) are listed below. 

1. The nature of the barrier in both grain boundary and c-axis junctions is unknown. 
In the better studied case of grain boundary junctions, we believe that the barrier is in 
the tunnel limit and for sufficiently large misorientation angles there are no pinholes in 
it. However, the typical barrier width and height don't seem to be well known. Therefore 
also the impact angle dependence of the barrier transmission is not known a priori. On 
the other hand, in order to explain the current-phase relation of 45° grain boundary 
junctions, we had to assume a weak impact angle dependence. It remains to be seen 
whether this agrees with the barrier properties determined independently, e.g. making 
use of the STM microscopy where the junction is viewed from above along a path crossing 
the grain boundary. With such a method, both the barrier height and width could be 
measurable. 

2. The role of various types of disorder in the junction region (interface roughness 
and faceting, spatially fluctuating barrier height and width, etc.) has to be studied sys- 
tematically. In particular, disorder is believed to reduce the ability of 45° grain boundary 
junctions to support mid-gap states [|39| and therefore diminishes the second harmonic 



as well. On the other hand, simple symmetry arguments predict a suppression of the 
first harmonic of I((p) by the surface roughness. Also an explicit calculation of the first 
harmonic in the presence of a finite barrier roughness supports this conclusion [ID|. A 
reliable answer to the question about which of the above two effects of disorder dominates 
is therefore crucial in order to decide whether disorder increases the ratio of the second 
and first harmonics of the current-phase relation \I 2 /Ii\ as suggested in Section 4, or it 
rather diminishes it. If the latter alternative is realized, then the experimental observa- 
tion of a large second harmonic provides even stronger argument in favour of anomalous 
quantum phase fluctuations in the (bulk) cuprates. 

3. Interface roughness and disorder in the barrier region may be also responsible for 
the weak impact angle dependence of the barrier transmission, which is implied by the 
large second harmonics observed in 45° grain boundary junctions. 

4. From the point of view of applications, especially the grain boundary Josephson 
junctions have attracted a lot of interest recently. In particular, junctions with sufficiently 
large second harmonics support doubly degenerate ground states (see Fig. [16]) and it has 
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Figure 16: Free energy F(ip) as a function of the phase difference across the weak link 
(taken from [||). The zero of energy has been set so that -F(O) = 0. (a) Asymmetric 
45° grain boundary. Top to bottom curves correspond to T = 40, 30, 20, 15, 10, and 4.2 
K, respectively, (b) Symmetric 45° grain boundary. Top to bottom curves correspond to 
T =20, 15, 11, 10, 5, and 1.6 K, respectively. 



been suggested |53| that this property might be exploited in the construction of a 'quiet 
qubit', i.e. of a two level system which couples only weakly to its environment and can 
be used as a basic element in a quantum computer. This is a fascinating proposal and 
a lot of efforts is being spent on its realization. Several nontrivial problems have to be 
solved before the final goal can be reached: on the technological side, a technique for a 
reproducible fabrication of well behaved submicron grain boundary junctions has to be 
developed. From the point of view of basic science, new routes to reducing the coupling 
of the junction to the environment are to be looked for, in order to make the junctions 
really quiet. 
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